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We consider a two- level quantum system (qubit) which is 
continuously measured by a detector and calculate the spec- 
tral density of the detector output. In the weakly coupled 
case the spectrum exhibits a moderate peak at the frequency 
of quantum oscillations and a Lorentzian-shape increase of 
the detector noise at low frequency. With increasing coupling 
the spectrum transforms into a single Lorentzian correspond- 
ing to random jumps between two states. We prove that the 
Bayesian formalism for the selective evolution of the density 
matrix gives the same spectrum as the conventional master 
equation approach, despite the significant difference in inter- 
pretation. The effects of the detector nonideality and the 
finite-temperature environment are also discussed. 



I. INTRODUCTION 

The long-standing and still controversial problem of 
quantum measurements!! is gradually becoming a prac- 
tical issue due to the development of experimental tech- 
niques capable_o£ measurements at the quantum border 
(see, e.g. Refs.BiD). The renewed interest in this subject 
is caused by the needs of quantum computing,cl since the 
measurement of an entangled and possibly evolving quan- 
tum state by a realistic detector in a realistic environment 
presents a nontrivial problem. _ 

Despite the experimental proofi of the violation of 
Bell's inequalityQ that rejects the idea of hidden vari- 
ables, the origin of the randomness of a quantum mea- 
surement result (the problem known as the "wavefunc- 
tion collapse") remains controversial. Important insight 
into this problem has been provided by the development 
of the theory of continuous quantum measurement, which 
generalizes the "orthodox" caseE2l of instantaneous mea- 
surement. There are tap different theoretical approaches. 
In the first approachO based on the-theory of interac- 
tion with a dissipative environment ,t3 the evolution of 
the density matrix of the measured system is averaged 
over a complete ensemble of measurements, thus leading 
to the deterministic equation. This is the best-known 
approach, at least in the solid-state community, and. aa 
can be called "conventional" . The other approachEj EZl 
(more developed in quantum optics) studies the random 
evolution of the quantum state during a single realiza- 
tion of the measurement, so that this evolution is con- 
ditioned on (selected by) the particular measurement re- 
sult. RecentlyB'ta the latter approach has been intro- 
duced into the context of solid-state physics using the 



new derivation based on the simple Bayesian analysis of 
probabilities. 

In the present paper this Bayesian formalism is applied 
to the calculation of the spectral density of the detector 
current when a two-level quantum system (qubit) is mea- 
sured continuously. As a particular example, we consider 
a double-quantum-dot occupied by one electron, the lo- 
cation p£adiich is measured by a quantum point contact 
nearbyEJ'Ea Another example of the setup is a single- 
Cooper-pak.pb.as, being measured by a single-electron 
transistor .tiEEa One more possible example is the con- 
tinuous measurement of two flux states of a SQUID by 
another inductively coupled SQUID.Ej 

We show that in the weak-coupling case when the 
quantum (Rabi) oscillations of the qubit state are only 
slightly perturbed by the detector, the corresponding 
peak in the spectral density of the detector current can 
be up, to 4 times higher than the noise pedestal (see also 
Ref.Ell As the coupling increases, the Quantum Zeno 
effected becomes significant leading to the Lorentzian 
shape of the spectral density centered at zero frequency. 

It is important to notice that there should be no dif- 
ference between the predictions of the conventional ap- 
proach and the approach of selective evolution unless the 
measurement result affects the system evolution (for ex- 
ample, via the feedback loop). We prove this equiva- 
lence explicitly for the detector spectral density (if there 
is no feedback). Despite the same final result, the in- 
terpretations are different: in the Bayesian formalism a 
significant contribution to the spectrum comes from the 
correlation between the detector noise and the system 
evolution, while this correlation is absent in the conven- 
tional approach. In the paper we also discuss the exten- 
sion of the Bayesian formalism to the case of additional 
weak interaction of the two-level system with a finite- 
temperature environment .E-j 



II. BAYESIAN FORMALISM 

For a two-level quantum system described by the 
Hamiltonian 

Wo = | (4* - 4 C 2) + H(c\c 2 + cjci) (1) 

(where e is the energy asymmetry and H is the tunnel- 
ing strength) the evolution of its density matrix pij due 
to continuous, measurement is given in the conventional 

chBE3e mi i 



approact 



by the equations 
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FIG. 1. Schematic of a qubit continuously measured by a 
detector with output signal I(t). 
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where the dephasing rate T due tc 
on the type of the detector used.l 

Eqs. (§)-(§) describe the averaged evolution. In con- 
trast, to analyze the single measurement process we need 
the selective (conditional) evolution of py which in the 
Bayesian formalism is described by equations^ 
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where I(t) is the detector current, Io = (Ii + /2V2, I\ 
and I2 are the average currents corresponding to two lo- 
calized states of the qubit, AI = I\ — (notice the 
different sign in the definition used in RefJIS), and Sq is 
the low-frequency spectral density of the detector shot 
noise (which is assumed to be constant in the frequency 
range of interest and to be practically independent of the 
qubit state) . The detector nonideality is described by the 
dephasing contribution 7 = T — (AI) 2 / 4So > due to in- 
teraction with "pure environment" (which does not affect 
the detector current). In particular, since T — U 
for symmetric quantum point contact (see Refs.1 
it represents an ideal detector,E3 r\ = 1, where r\ = 1— 7/T 
is the ideality factor. In contrast, the single-electron tran- 
sistor in a typical operaijpouDoint is a significantly non- 
ideal detector, 77 <C l.EfBEj (Actually, for the single^, 



n 

electron transistor Eq. (^|) can be further improved,c 
however, it does not make a difference for the purposes 
of the present paper.) The SQUID can be aijt-jdeal de- 
tector when its sensitivity is quantum-limited.E3'E!J Since 
the typical output signal from the SQUID is the voltage 
(not the current), this requires a minor modification of 
the formalism; so in this paper we do not consider the 
SQUID case explicitly even though all the results can be 
easily translated into SQUID language. 

Eqs. (@)-@ allow us to calculate the evolution of py 
if the detector output I(t) is known from the experi- 
ment. In ordeiy-ko simulate the measurement we need 
the replacement!!^! 



I(t)-Io = AI(pn-p 22 )/2 + £(*), 



(6) 



where the random process £,(t) has zero average and 
"white" spectral density S$ = So- 
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FIG. 2. Spectral density Si(uj) of the detector current for 
weak coupling (a — 0.1) of an ideal detector (rj = 1) with 
a two-level system (e/H — 0, 1, and 2). The results of the 
Monte-Carlo calculations using Bayesian formalism are shown 
by solid lines, while the dashed lines (almost coinciding with 
solid lines) are calculated using the master equation approach. 



Notice that in order to consider the detector as a 
device with classical output signal, Eqs. (|^)~(||) essen- 
tially use the Markov approximation and assume that 
the typical frequency of the internal detector evolution 
(on the order of So/e 2 ~ lo/e) is much higher than 
the typical frequency max(r2,T) of py evolution (here 
Q = V 4H 2 + e 2 /h is the frequency of unperturbed quan- 
tum oscillations). In particular, this condition requires 
the detector to be "weakly responding", |A/| <§; I , that 
allows us to use the linear response approximation. 

Averaging of Eqs. (|J)-(H) over all possible measure- 
ment results (i.e. over random contribution £(t)) reduces 
them to Eqs. (§)-(||). Notice that the stochastic equa- 
tions are written in Stratonovich form which preserves 
the usual calculus rules, while averaging would be easier 
in Ito formH3 



III. WEAK COUPLING 

Using Eqs, (j^)-(||) and the Monte-Carlo method (sim- 
ilar to Ref.ta) we can calculate in a straightforward way 
the spectral density Si(uj) of the detector current I(t). 
Solid lines in Fig. || show the results of such calcula- 
tions for the ideal detector, rj = 1, and weak coupling 
between the qubit and the detector, a — 0.1, where 
a = h(AI) 2 /8S0H (a is JL times less than the param- 
eter C introduced in Ref.Ej). One can see that in the 
symmetric case, e = 0, the peak at the frequency of quan- 
tum oscillations is 4 times higher than the noise pedestal, 
Si (fl) = 5So while the peak width is determined by the 
coupling strength a (see Fig. |^ below) . In the asymmetric 
case, £ 7^ 0, the peak height decreases (Fig. ||), while the 
additional Lorentzian-shape increase of Si (to) appears at 
low frequencies. The origin of this low-frequency feature 
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FIG. 3. A particular realization of the evoltion of pu(t) 
due to continuous measurement for e/H — 1, a = 0.1 and 
7] = 1. Notice the fluctuation of both the phase and the 
asymmetry of oscillations. 

is the slow fluctuation of the asymmetry of pu oscilla- 
tions (Fig. H). In case e = the amplitude of pu oscilla- 
tions is maximal (see thick line in Fig. [|a), hence there is 
no such asymmetry and the low-frequency feature is ab- 
sent, while the spectral peak at the frequency of quantum 
oscillation is maximally high. 

In order to understand the factor 4 for the maximum 
peak height, let us consider the case a -C 1, e = 0, and 
rj = 1. Then the selective evolution can be written as the 
quantum oscillations with slowly fluctuating phase pit): 

z{t) = p n {t)-p 22 (t)= cos <f>(t), cf> = nt + tp(t), (7) 
p 12 = {i/2)sm<f>(t) (8) 

(the state becomes pureEl after an initial transient period 
since r\ — 1, while the real part of p 12 decays because of 
e = 0). From Eqs. (|) and (§) we obtain the random 
phase dynamics: 



ip = -sin0^[cos^^ + C(t)]- 



(9) 



Since (AI) 2 /2So -C ^, we can neglect the first 
term in the square brackets and average the sec- 
ond contribution over (f> that leads to the white 
spectrum of <p: S^iuj) = (AI) 2 /2S . Then the 
correlation function K z (t) = (z(t)z(t + r)) can 
be calculated as K z (t) = cos(17r)(cos A<p(r))/2 = 
cos(Qt) exp[— (AI) 2 t/8Sq]/2 and the spectral density 
S z (lo) = 2 J_ K z (t) exp(iu)T)d,T has a peak in the vicin- 
ity of the oscillation frequency, \u> — (l\ <C f2: 



S z (uj) = 
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The detector current is given by Eq. (|6|), so its spectral 
density consists of three contributions: 



Si(w) = So 



AI 2 



-S z (u) 



AI 



SfzM + S^M], (11) 




FIG. 4. A particular realization of pu evolution (thick line) 
and the corresponding detector output I(t) (thin solid, dot- 
ted, and dashed lines) averaged using rectangular windows 
with different time constants r a . (a) - weak coupling case, 
a = 0.1, (b) - strong coupling case, a — 5. 



where the last contribution originates from the cor- 
relation between evolution and the detector noise 
£(i). To calculate the correlation function K^ z (t) = 
(£(t)z(t + t)) for t > 0, we need to take into ac- 
count the phase shift dtp = — sin (fiAIS^^fydt dur- 
ing even infinitesimally small time duration dt, since 
the amplitude of the stochastic function grows 
with the timescale decrease, ^(t) 2 dt = So/2. Us- 
ing trigonometric formulas and linear expansion in 
dp we obtain (£(t) cos[(p(t) + d(f> + fir + Ap(r)]) = 
A.ISo 1 {£(t) 2 dt) (sin <j){t) sin[(f)(t)+ClT]) (cos Ap(r)) and fi- 
nally K (z (t) = A/cos(^T)exp[-(A/) 2 T/8S'o]/4. After 
the Fourier transformation one finds that the correlation 
between £(i) and z(t) brings exactly the same contribu- 
tion to the detector spectral density (see Eq. (|TT]) ) as the 
term due to z(t) evolution, so that 



Si{bj) = S + 



4S 



1 + [85 (w - fi)/(AJ) 



212 



(12) 



Thus, the peak corresponding to quantum oscillations is 
4 times higher than the noise background, while its full 
width at half height is equal to (AI) 2 /4Sq =j-Q^ (the 
same peak width has been calculated in Ref£3). The 
integral under the peak, 



27T 4 



(13) 
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has an obvious relation to the average square of the detec- 
tor current variation due to oscillations in the measured 
system. Notice, however, that this integral is twice as 
large as one would expect from the classical harmonic 
signal, since one half of the peak height comes from non- 
classical correlation between the qubit evolution and the 
detector noise. Classically, Eq. would be easily un- 
derstood if the signal was not harmonic but rectangular- 
like, which is obviously not the case. Actually, the detec- 
tor current shows neither clear harmonic nor rectangular 
signal distinguishable from the intrinsic noise contribu- 
tion. Figure ^a shows the simulation of pn evolution 
(thick line) together with the detector current I(t). Since 
I(t) contains white noise, it necessarily requires some av- 
eraging. Thin solid, dotted, and dashed lines show the 
detector current averaged with different time constants 
r a : T a fl/2n — 0.3, 1, and 3, respectively. For weak aver- 
aging the signal is too noisy, while for strong averaging 
individual oscillation periods cannot be resolved either, 
so quantum oscillations can never be observed directly by 
a continuous measurement (although they can be calcu- 
lated using Eqs. (f|)-(|5])). This unobservability is revealed 
in the relatively low peak height of the spectral density 
of the detector current. 



IV. ARBITRARY COUPLING 
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The situation changes as the coupling between the de- 
tector and qubit increases, a > 1. The strong influ- 
ence of measurement destroys quantum oscillations, and 
the Quantum Zeno effectH develops, so that for a 3> 1 
the qubit performs random jumps between two localized 
states (see Fig. In this case the properly averaged 
detector current follows pretty well the evolution of the 
qubit (however, the unsuccessful tunneling "attempts" 
still cannot be directly resolved), and the spectral den- 
sity of I(t) can be_calculated using the classical theory 
of telegraph noisecfl leading to the Lorentzian shape of 
Si (uj) . Figure |^a shows the gradual transformation of the 
spectral density with the increase of the coupling a for 
a symmetric qubit, e — 0, and an ideal detector, rj = 1. 
The results for an asymmetric qubit, e/H = 1, are shown 
in Fig. |b. 

The curves in Fig. ||as well as the dashed curves in Fig. 
^ are calculated using the conventional master equation 
approach which gives the same results for the detector 
spectral density as the Bayesian formalism (we will prove 
this later). In the conventional approach we should as- 
sume no correlation between the detector noise and the 
qubit evolution (the last term in Eq. ( |ll| ) is absent) while 
the correlation function (r) should be calculated con- 
sidering z(t) not as an ordinary function but as an oper- 
ator. Then the calculation of (z(t)z(t + r)) can be essen- 
tially interpreted as follows. The first (in time) operator 
z(t) collapses the qubit into one of two eigenstates which 
correspond to localized states, then during time r the 



FIG. 5. The detector current spectral density Si(uj) for 
rj — 1 and different coupling a with (a) symmetric (e = 0) 
and (b) asymmetric (e/H = 1) qubit. 

qubit performs the evolution described by conventional 
Eqs. (g)-(||), and finally the second operator z(t+r) gives 
the probability for the qubit to be measured in one of two 
localized states. j-XOf course, this procedure can be done 
purely formally^ without any interpretation.) Notice 
that there is complete symmetry between states "1" and 
"2" even for e ^ (in particular, in the stationary state 
Pn = P22 = 1/2), so the evolution after the first col- 
lapse can be started from any localized state leading to 
the same contribution to the correlation function. In this 
way we obviously get Kz(t) = Ph(t) — Pti(t) where pa 
is the solution of Eqs. (§)-(||) with the initial conditions 
,011(0) = land pi 2 (0) = 0. 

For the symmetric qubit, e = 0, these equations can 
be easily solved analytically, and finally we obtain: 

Q 2 (AI) 2 T 

s ^> = 5 ° V-V+r^ ' (14) 

where T — rj^ 1 (AI) 2 /ASq — cm - 1 SI. This equation obvi- 
ously transforms into Eq. ( |12| ) for rj — 1 and a < 1. 
Notice that for weak coupling with a nonideal detec- 
tor, rj < 1 and arf 1 <C 1, the peak height of Si(uj) is 
equal to ArjSo, while the linewidth arj^ 1 ^ of the peak 
is r; _1 times wider than for the ideal detector. As the 
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coupling increases, the linewidth grows and the oscilla- 
tion frequency decreases^ lo osc — f2[l — (a/277) 2 ) 1 / 2 ]. 
The transition into the overdamped regime occurs at 
otr\~ x > 2 while the peak- like feature disappears at 
y/2. For arj 1 > 2 the spectral density con- 
sists of two Lorentzians [w 1)2 = T/2 =F (T 2 /4 - ft 2 ) 1 / 2 ] 
centered at zero frequency, with the negative sign and 
the smaller amplitude Ai of the second Lorentzian which 
has higher cutoff frequency, A%jA\ — —<jj\/u)2- In the 
case ai]^ 1 ^> 1, which corresponds to the well developed 
Quantum Zeno effect, Si(uj) — Sq has purely Lorentzian 
shape (A/) 2 wi/(lj 2 + uj 2 ) with w\ = fi 2 /F = flrj/a. 

For the asymmetric qubit, e ^ 0, the spectral density 
can in principle also be calculated analytically but the 
expressions would be too lengthy and it is simpler to 
use numerical solution of Eqs. (||)-(||). The analytical 
formula for the weak coupling limit is 



vS e 2 /H 2 



l + (uh 2 n 2 /AH 2 T) 2 
4t 1 S (1 + s 2 /2H 2 )- 1 



i + [(uj -n)/r(i-2H 2 /n 2 n 



2o2M2 



(15) 



The spectral peak and the low-frequency Lorentzian be- 
come wider with the coupling increase since F = a?7 -1 f2, 
and for \e/H\ < l/v2 the overdamped regime starts 
from T = Ti where T 2 12 = (fl 2 /2a)[b =F {b 2 - 4a) 1 / 2 ], 
b ee 1/4 - 27a 2 /4 + 9a/2, a = e 2 /(AH 2 + e 2 ). At 
r > T2 the dynamics formally returns to the under- 
damped regime, however, the peak linewidth is much 
larger than the frequency and so Si(u>) is monotonous. 
For \e/H\ > l/v2 the overdamped regime does not oc- 
cur. In both cases in the limit of large V the spectral 
density has almost Lorentzian shape with the cutoff fre- 
quency ui = AH 2 /fi 2 T. 

One can check that the spectral densities given by Eqs. 
( |l4| ) and (|l5|) satisfy the integral condition (|l3|), which 
remains valid for arbitrary parameters a, e/H, and 77, 
because of the equation Ki(+0) = (AI/2) 2 . 



V. EQUIVALENCE OF TWO APPROACHES 

Comparing two derivations of Si(u>) in the case a <C 1, 
i] = 1, and e = 0, we see that K%{t) is twice larger than 
K z (t) because in the conventional approach the evolu- 
tion always starts from the localized state while in the 
Bayesian approach it starts from arbitrary phase of quan- 
tum oscillations. This difference exactly compensates for 
the absence of the last term in Eq. ( ]Tl| ) in the conven- 
tional approach. 

Let us prove explicitly that the two approaches give 
the same result for Si(u>) in a general case. In order to 
calculate K^ z (t) for t > using the Bayesian formalism, 
let us first average the product £(to)z(to+r) over random 
during time period to < t < to + t, fixing the same 
conditions at t = to. Then we can use conventional Eqs. 



d)-@ [regarded as Eqs. averaged over random 

with the initial condition pij(to + 0) = Pij(to) +dpij 
where 

dz = AIS - 1 [l-z 2 {t )]£(t )dt, (16) 
d Pl2 = -AIS - 1 z(to)pi 2 (to)Z(to)dt (17) 

(for simplicity we will refer to z = pn — P22 as a 
component of p%j). Since the sign of £(io) is arbitrary 
and averaged evolution equations are linear, we need 
only the fluctuating contribution to pij (to + 0) and, 
hence, can formally assume that the evolution starts from 
Pijito + 0) = dpij (notice that we can forget the condi- 
tion pn + P22 — 1 and use only z and p\2)- Using the 
relation ^{to) 2 dt = ^0/2 and the evolution linearity, we 
can formally write 

K iz (r) = (AI/2)(z(t + r)), (18) 

where pij satisfies Eqs. (H)-(||) with z(to) = 1 — z(to) 2 
and pi2(to) = —z(to)pi2(to), and the averaging over 
the initial conditions at t = to should still be done 
later. Before that let us do a similar formal trick for 
K z (t) representing it as (z(to + r)) where the evolution 
starts from S(t ) = z(t Q ) 2 and pi 2 (to) = z(t )pi2(to)- 
Now combining two terms in the detector current cor- 
relation function Kj(t) = (I(t )I(t + r)) - (I) 2 = 
(AI/2) 2 K z {t) + (AI/2)K (x (t) (here r > 0), we see that 
it can be written as (AI /2) 2 (z(to + r)) where z(to) = 1 
and pi2(io) — 0. Thus we have exactly arrived at the 
expression of the conventional approach, in which the 
evolution always starts from the localized state, regard- 
less of the actual quantum state at t = to- This proof 
is obviously valid for arbitrary a, 77, and e/H. Despite 
the same result in two approaches, the interpretations 
are quite different since the Bayesian approach allows 
us to follow the qubit evolution during the measurement 
process, while the conventional approach gives only the 
average characteristics. 



VI. FINITE-TEMPERATURE ENVIRONMENT 

In this section we will discuss how to introduce the 
finite-temperature environment into Eqs. (Q)~(|(]) of the 
Bayesian formalism. Notice that so far there has been 
complete symmetry between the states "1" and "2" even 
for finite energy difference e, while the finite tempera- 
ture effects would be expected to lead to different aver- 
age populations of these states. Such symmetry requires 
an implicit assumption that the typical energy in the de- 
tector (voltage or temperature) is much larger than the 
energies involved in the qubit dynamics. So, the absence 
of temperature in the formalism does not mean that it is 
zero or very large, it just means that the temperature ef- 
fects are not important. Now let us assume that besides 
the detector, the qubit is coupled to an additional finite- 
temperature environment, which creates an asymmetry 
between states "1" and "2" when e 7^ 0. 
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While the case of finite coupling of a two-level system 
with an environment presents a really difficult problem,!! 2 ] 
the weak coupling limit-can be treated in a simple way. 
In the standard methods it is described by the equations 



P++ = -7l0++ -Pst), 

p + - = ifip+_ - 72P+-, 



= 1, 



(19) 
(20) 



which are written in the diagonal basis ( "+" corresponds 
to the ground state) . The temperature T determines the 
stationary occupation p st = [1 + exp(— Ml/T)]^ 1 of the 
ground state, and the relaxation rates obey inequalityE3 

71/2 < 72 < n. 

If the coupling of the qubit with the detector is also 
weak, oa]~ x <C 1, the evolution due to extra finite- 
temperature environment can be simply added to the 
evolution due to measurement. For this purpose Eqs. 
(pj|)-(p0|) should be translated into the basis of localized 
states, so the terms 



-(A 2 7l + B 2 l2 ){ Pll - 1/2) - 7l A(l/2 - Pst ) 

+-4-B(7! - 7 2 )Repi2, 
where A = e/hfl, B = 2H/HQ, 

should be added into Eq. 



(21) 



for pn and the terms 

-{A 2 l2 + B 2 7 i) Repia + AB(p u - p 22 )("fi - 72)/2 
+7lB(l/2-p st )-j 72 Imp 12 (22) 

should be added into Eq. @ for p i2 . The same terms 
should obviously be added into Eqs. @-(||) for the con- 
ventional approach. (Of course, this generalization is 
purely phenomcnological and is limited to the weak cou- 
pling regime, so the effect of Eqs. (^l|)-(p2[) can be consid- 
ered only at the timescale longer than oscillation period.) 

In the generalized case it is still possible to prove 
that the results of the Bayesian formalism for the detec- 
tor current spectral density Si(u>) exactly coincide with 
the results of the conventional approach. The essen- 
tial difference from the proof above is nonzero station- 
ary solution (z st , pi2,st) °f modified Eqs. (||)-(§) when 
Pst 7^ 1/2. It is convenient to consider homogeneous evo- 
lution equations (with p st — 1/2) simply shifting z(t) 
and pi2{t) by the stationary values. Using the same 
idea as in the proof above we can show that in the 
Bayesian approach Kj(t) for r > can be written as 
(A//2) 2 (z(tn + r)) where pij satisfies homogeneous mod- 



ified Eqs. (0)-© with z(t ) = 1 - 2z st z(t ) + zj t and 
pu(to) = -z(t )pi2, s t - z st pi2(t ) + z st pi2,at- After the 
averaging over initial states these initial conditions can 
be replaced with z(to) = 1 — z 2 st and p~\i{to) = —z s tPxi,st- 
Now let us show that we get the same Ki(j) in the 
conventional approach. With the probability (1 + z st )/2 
the first operator z(to) localizes the qubit into the state 
"1". Then the initial state for the homogeneous equa- 
tions is z(to) = 1 - z st , pxiito) = -p\2,st- With the 
probability (1 — z s t)/2 the evolution starts from the state 
"2", i.e. z(t ) = -1 - z st , Pi 2 (*o) = ~Pi2,st- Adding two 



contributions with opposite signs we see again that for 
t > 0, Kj(t) = {AI/2) 2 z(t + t) where z can be found 
as a solution of Eqs. @-(||) modified by Eqs. (21)— (^S 



without inhomogeneous terms, with the initial condition 
z(t ) = 1 - z 2 st and pi2(h) = ~z st p 12 , s t- Thus, the cor- 
relation function Ki(t) and, hence, the spectral density 
Si{u>) coincide in the two approaches. 

It is technically simpler to consider the averaged evo- 
lution in the diagonal basis rather than in the basis of 
localized states [for this purpose we need to translate the 
term —Tpi 2 from Eq. (|^) into the diagonal basis and add 
it into Eqs. (fl9|)-(|2C|)]. So, to calculate the correlation 
function Kg(r) analytically, we start the evolution from 
one of the localized states, then consider the averaged 
evolution in the diagonal basis (neglecting the rapidly 
oscillating terms due to measurement), and make the sec- 
ond projection onto localized states at t — r. Finally we 
obtain the resulto 



Si(u>) = Sq + 



(A/) 2 



where 



Zst 



W t 

2(AI) 2 H 2 1 

w h 2 w 1 + [(w - n)/w } 2 



_£ 1 m 

i + 4ip r / 7l a 2 ^ tanh( 2 r } ' 



{w/W t f 



W t = 7 i 



ATH 1 

h 2 n 2 ' 
r 

Wo = 7 2 + -(l + 



n 2 n 2 



(23) 

(24) 

(25) 
(26) 



Let us emphasize that the effect of finite-temperature 
environment is not generally equivalent to the nonideality 
of the detector described by finite 7 in Eq. (^|). As an 
example, in the case of extra environment the right hand 
part of Eq. (|l^) should be multiplied by the factor 1 — z 2 t 
which disappears (z st = 0) only if T = 00 or e = 0. 

Comparing Eqs. (|23| ) and ( |l5| ) one can see that within 
the accuracy of weak coupling approximation the change 
of Si(uj) due to extra environment can be reduced to the 
detector nonideality, 77 < 1, in two cases. If \e/H\ <C 1, 
then Eqs. (|2^) and ([l5]) coincide at arbitrary tempera- 
ture T for 77 = (1 + 272 /ail)" 1 . For asymmetric qubit, 
\s/H\ > 1, the equivalence is possible only at high tem- 
peratures, T ^> HQ, and requires conditions 



72 =7i(l + £ 2 /2# 2 )/2, 

rT 1 = 1 + (1 + e 2 /AH 2 f ,2 ll /an. 



(27) 
(28) 



Figure |6J shows the numerically calculated spectral den- 
sity Si(u) of the detector current for a nonideal detector, 
i] = 0.5 (dashed lines) and for an ideal detector but extra 
coupling of the qubit to the environment at temperature 
T = H (solid lines). The rates 71 and 72 are chosen ac- 
cording to Eqs. ( p7| ) and (|2"s|). For the symmetric qubit, 
e = 0, the results of two models practically coincide. In 
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FIG. 6. Spectral density Si(u>) for the nonideal detec- 
tor (dashed lines) and for the case of the ideal detector and 
weak extra coupling with the finite-temperature environment, 
T = H (solid lines). 



oscillations. This makes it possible to tune the phase us- 
ing the feedback control of the qubit parameters. If the 
real-time calculations using Eqs. (@)-(§) and fast feed- 
back loop were available in an experiment (the typical 
bandwidth should be larger than T) then the random 
diffusion of the oscillation phase could be eliminated and 
the qubit could "stay fresh" for a very long time. The 
suppression of qubit dephasing using the feedback control 
of the tunneling strength H has been confirmed by the 
Monte-Carlo simulations. The elimination of the phase 
diffusion gives rise to the (5-function peak in the detec- 
tor spectral density Si(u>) at the frequency Q. A de- 
tailed analysis of this situation is beyond the scope of 
the present paper. 

The author thanks D. V. Averin, J. E. Lukens, and K. 
K. Likharev for fruitful discussions. The work was partly 
supported by AFOSR. 



contrast, the solid and dashed lines for e = 2H signif- 
icantly differ from each other at low frequencies, while 
the spectral peak at u) ~ £1 is fitted quite well. 



VII. CONCLUSION 

Using both Baycsian and conventional approaches, we 
have calculated the spectral density Si{u>) of the detector 
current when a two-level quantum system (qubit) is mea- 
sured continuously. Depending on the coupling strength, 
there is a gradual transition from quantum oscillations 
to quantum jumps. This results in a transition from 
the peak-like spectral density to the Lorentzian shape of 
Si(u>). The maximal height of the peak at the frequency 
of quantum oscillation is shewn to be 4 times the shot 
noise pedestal (see also Ref.cj). 

In the simple case of weak coupling between a symmet- 
ric qubit and an ideal detector, the height of the spectral 
peak is twice as high as the classical result for a har- 
monic signal. In the Bayesian approach this is explained 
by the significant correlation between the detector noise 
and the further evolution of the measured system due to 
quantum back-action. In contrast, in the conventional 
approach this fact is the consequence of discrete eigen- 
values of z operator, which corresponds to the magnitude 
measured by the detector. (In other words, this opera- 
tor "collapses" the system into one of two eigenstates, 
and that is why the averaged product of two operators 
is twice as large as that for a classical harmonic signal.) 
So, even though the results for Si(u>) coincide in two 
approaches, the interpretations are significantly different 
since the "abrupt" collapse is replaced in the Bayesian ap- 
proach by the "continuous" collapse related to the noisy 
detector output. 

It is important to notice that the Bayesian formalism 
allows us to monitor continuously the phase of quantum 
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